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Abstract 

One-dimensional, boundary-driven lattice gases with local interactions are stud- 
ied in the weakly interacting limit. The density profiles and the correlation func- 
tions are calculated to first order in the interaction strength for zero-range and 
short-range processes differing only in the specifics of the detailed-balance dy- 
namics. Furthermore, the effective free-energy (large-deviation function) and the 
integrated current distribution are also found to this order. From the former, we 
find that the boundary drive generates long-range correlations only for the short- 
range dynamics while the latter provides support to an additivity principle recently 
proposed by Bodineau and Derrida. 
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1 Introduction 



One-dimensional boundary driven lattice gases are simple model systems which 
allow detailed studies of nonequilibrium steady-states. They are relevant in the 
sense that steady states in experiments are often produced by choosing appropri- 
ate boundary conditions (heating a horizontal layer of liquid from below being 
the most widely quoted example) and, furthermore, these models appear to cap- 
ture some important consequences of the nonequilibrium drive such as e.g. the 
generation of long-range correlations. 

Analytical approaches to driven lattice gases proceed along two distinct paths. 
On one hand, exploiting the matrix product method, Derrida and coworkers [HIH 
El |4j |6) have come up with a series of exact results for the exclusion process 
(hard-core particles undergoing symmetric or biased diffusive motion) in one di- 
mension . On the other hand, in a parallel series of papers, Bertini and coworkers 
00 solved some of the same problems employing the formalism of fluctuat- 
ing hydrodynamics. The results obtained for the probability functional of a given 
density profile or for the integrated current distribution are important since they 
provide new insights into the properties of nonequilibrium steady states. Thus it 
would be highly desirable to establish how general or universal these results are. 
In order to investigate this generality issue, in this paper, we go beyond the speci- 
ficity of the exclusion process, and study the effects of varying the interactions 
between the particles as well as varying the microscopic dynamics in boundary 
driven lattice gases. 

Our main "technical" results are the derivation of the weak interaction limits of 
both the nonequilibrium free energy and the distribution of the integrated particle 
current. From the explicit form of the free energy, we can then deduce the density 
profile and the correlations (effective interactions) for various interactions (pair or 
triplet) and various dynamics [zero-range [ 1 1 and short-range (misanthropic ifTTI ) 
processes]. 

Our main "physical" findings are that the density profiles are nonlinear and 
they depend on both the interactions and the dynamics. Furthermore, the details 
of the dynamics are found to change the the correlations from short- to long-range, 
and they are shown to be able to change the sign of the effective interactions. As 
to the distributions of the integrated particle current, we verify that, for all the 
interaction- and dynamics combinations we studied, they are in agreement with 
the recently proposed addititvity principle which allows to construct this distribu- 
tion function from the knowledge of its first two moments lfT2l. 

We start by introducing the models of interacting lattice gases and deriving dy- 
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namical rules satisfying detailed balance in the absence of boundary drive (Sec©. 
Next we describe on the example of free particles how the models can be formu- 
lated in terms of field theories (Sec|3]). Zero-range and short-range (misanthropic) 
processes are worked out in Sections |4] and respectively. There we provide 
explicit expressions for the effective free energies. Finally, Sec|^lis devoted to a 
separate treatment of the integrated current distribution. 

2 The model 

2.1 Lattice gas with onsite interaction and boundary drive 

A d — 1 dimensional lattice gas is considered with hopping dynamics in the 
bulk and with particle injection and removal at the ends of the chain. The state 
of the system n = {n , n 1 , n L } is specified by the number of particles rii = 
0, 1, oo at lattice sites i (i = 0, 1, L), and the interactions are assumed to be 
local 

L 

E (n) = Y J h i (1) 

i=0 

where hi = h(rii) is the energy of interactions among particles at the same site. 
We shall treat the simplest case of pair interactions, h{n) = en(n — 1). However 
triplet interactions, h{n) = en(n — 1) [n — 2), will also be considered occasionally 
either for demonstrating the effects of competing interactions 

hi = h{rii) = e (±nj(nj - 1) + An^rij - l)(nj - 2)) . (2) 

or for probing the universality of our results with respect to varying the type of 
interaction. For A = 0, we restrict our study to the + sign in the r.h.s. of © in 
order to have stability against collapse of all particles on a single site, while for 
A > both signs will be considered. 

The dynamics of the system is described in terms of a master equation for the 
time evolution of the probability P(n, t) of a given particle configuration 

d t P(n, t) = C D P(n, t) + C BC P(n, t) (3) 

where the first and second terms on the right-hand side represent the nearest- 
neighbor hopping in the bulk and the injection-removal terms at the boundaries. 
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The bulk term is given through the rate of hopping w i ^ i ± 1 (n) from site i to i ± 1 
in a state n as 

L-1 L 

C D P(n,t) = -J2^i+i(n)P{n,t) -Y,^^~i(n)P{n,t) (4) 

j=0 i=l 
L-1 L 

+ iu i _ >i+ i(n i+ i_ >i )P(n i+ i_ >i , t) + ^2 w i ^ i - 1 {n i ^ i )P{n i _ 1 ^ t) 

j=0 i=l 

where state is obtained from n by moving a particle from i to j. The injection 
and removal of particles at the boundaries % — and 2 = L are described by the 
terms 

C BC P(n,t) = - [iu+(n) + u>+(n) +^(n) + u£(n)] P(n,t) (5) 
+ wJ(no-) P(n -,t) + u>o~(n +) P(n +,t) 
+ wl (n L - ) P{n L - ,t)+wl (n L+ ) P(n L + , i) 

where u^n), Wq(u), w~^(n) and w^(n) are the rates of adding (+) or removing 
(— ) a particle at site or L in the state n and, furthermore, the state n i+ (_) is 
obtained from n by adding (removing) a particle at site i. 

2.2 Choice of dynamics 
2.2.1 Equilibrium distribution 

In order to specify the dynamics, we shall assume that if the injection and removal 
rates are such that there is not net flux (particle or energy current) through the 
system then equilibrium is reached. Due to the onsite nature of the interaction, the 
equilibrium distribution factorizes, it becomes the grand-canonical distribution 

1 (^ n iQ~ h i/ T 

P sq (n) = [[p eq (ni), PeqK) = ^ j (6) 

i 

where T and ( are characterizing the temperature and fugacity of the outside heat 
and particle bath, and S is a normalization factor. In the following, we shall absorb 
T in the interaction strength e and the high temperature limit studied below will 
simply mean that e is small. 
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2.2.2 Hopping rates 



The bulk hopping rates may now be defined by making the natural (though not 
necessary) assumption that they satisfy detailed balance with P eq , and that the 
detailed balance remains valid in the presence of the boundary drive, as well. This 
assumption yields the following condition for the rate of hopping to the right 

Wi^ i+1 (n) n . e -[E(n^ l+1 )-E(n)]/2 



W i+ l_rf(7^ + i) {U i+ i + l) e -[E(n)-E(n^ +1 )]/2 

„ ~h(ni)—h(rii—i) 

(8) 



(n i+ i + i) e ^(ni+i+i)-^K+i) 



and similar condition with i + 1 replaced by i — 1 applies for the rate of hopping 
to the left. Clearly, the hopping rates are not uniquely determined by detailed 
balance. The remaining arbitrariness is not relevant for equilibrium but, in the 
presence of a drive, the steady state properties do depend on details of the dynam- 
ics. In order to have some idea about the importance of various choices, we shall 
consider two significantly different rates. 

First, the rates will be chosen so that they depend only on the energy change 
at the site from which the hopping originates. As can be seen from ©, a choice 
satisfying this condition is as follows 

w\Z] +1 (n) = Dn^ n ^ h ^-^ (9) 

where D is the "diffusion coefficient" setting the timescale. Dynamics defined by 
the above rate is an example of the so called zero range processes lfTUl[P31 where, 
in general, the rate is an arbitrary function of n«. Thus, the first half of our study 
can be viewed as an investigation of zero range processes in the presence of a 
boundary drive. 

For a second choice, a more standard hopping dynamics will be used with rates 
depending not only on the energy change at the starting position of the hopping 
particle but on the total energy change due to the hopping. This choice follows 
from the eq.© version of detailed balance 

W {ST) - +l {n) = D ni Q-lH^- 1 )+Hni+i)-Hn i )-h{n i+1 )]/2 

where a subscript in signifies the name short range process ifTTl we shall 
use in order to distinguish the resulting dynamics from the zero range process. 
An important feature of the short range process is that tuning the interactions 
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(0 < e < oo, A = 0) so that one extrapolates between the exactly solvable nonin- 
teracting limit (e = 0) and the much investigated and also exactly solvable EE! 
boundary driven symmetric exclusion process (e — > oo). 



2.2.3 Boundary drive (injection and removal rates) 

Nonequilibrium drive is introduced into the model through injection and removal 
of particles at the ends (j = and L) of the chain. The rates can be set again by 
referring to detailed balance condition with the baths attached to the ends of the 
chain. The baths are assumed to contain noninteracting particles with the chemical 
potentials set to produce the appropriate injection and removal rates. 

For zero range process, the removal can be imagined as diffusion into the 
baths. The rates are independent of the properties of the baths and given by © 
with D replaced by new (externally controlled) parameters j3 and 7, yielding 

ioJ zr) ~(n) = inQQ Kn a )-h{n -i) . w ^~(n) = /3n L e h{riL) - h{nL - 1) . (11) 

Within the picture of zero range process, the injections (moving a particle from 
the bath to the end sites) are independent of state of the end sites, and thus they 
are constants a and 5 determined by the properties of the baths 

w^ )+ {n) = a ; w f )+ (n) = 5. (12) 

The boundary rates for the short range process are determined along similar 
lines. The rate of removal follows from and, since the particles in the baths are 
not interacting, they are given by 

4 sr) -(n) = inQQ H^Kn Q ~l)]/2 . w («)-(tf) = p nLfS M"r.)-*{nr.-W2. (1 3 ) 

In order to satisfy the detailed balance Q, the injection rates must also depend on 
the energy change due to the adding a particle and, using again that the particles 
in the bath do not interact, we find 

W ( ° T)+ (n) = ae Mn )-h(n +l)]/2 . w W+(fy = 6e [h(n L )-h(n L+ l)]/2 _ (U) 

We have now finished the description of the two models we shall be concerned 
with in this paper. Eqs. d9ll 11121) define the boundary driven zero range process 
while the rates given by eqs. (I10I13I14I) define the corresponding short range 
process. As one can easily verify, the case of 

-='* (15) 

7 P 
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corresponds to equilibrium with the baths (no drive) and the steady state distri- 
bution is the equilibrium distribution given by eq.©. Nonequilibrium drive is 
present when the above equality is violated. 



3 Preparing for the perturbation expansion: Non- 
interacting particles 

3.1 From the master equation to path integrals 

We shall now perform the mapping of the master equation for the probability that 
the system is in the occupation numbers configuration n at time t, namely, P(n; t), 
to a field theory, where approximation techniques inherited from our experience 
with equilibrium systems can be transferred without any difficulty. 

We only briefly sketch the procedure (see ifbUl for a recent review). We build 
up the state vector \^(t)) = J2f{P(n,t)\n) and rewrite the master equation for 
P in terms of and evolution equation for \^(t)). This equation can be cast in the 
form 

d|¥(t)> 



dt 



Co\*{t)) (16) 



where C is an evolution operator acting on the microstates \n) indexed by the 
local particle numbers. Not suprisingly, the evolution operator C is conveniently 
expressed in terms of the creation and annihilation operators a\, a* related to each 
local occupation number n». For the rates specified in the previous subsection, but 
at s = 0, it takes the following form 



(a] - a]) a* 



-A.=J>E £ 

i j nn of i 

+ a(l - aj) + 8(1 - 4) (17) 
+ 7(aJ - l)o + P(a\ - l)a L 

One can of course recognize that at e = one is formally dealing with a free boson 
hamiltonian. The next step consists in converting the computation of expectation 
values of physical observables into the evaluation of a path-integral over two fields 
ai(t) and cii(t). This step leads to averaging a physical observable A({ni}) in the 
following way: 

(A(n)} = J Va^ A(a(t)) e' 3 ^ (18) 
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where 



S [a,a\ = -y^gj(t) + y dr (did T ai + D ^(dj - %)q») 



+ a(l - a ) + 7(a - l)a 
+ 5(1 -a L ) +(3(a L - l)a L 



(19) 



Note that we have moved the initial condition to t = — oo, in order to sit directly 
in the steady-state (initial condition independent). The prescription to obtain A 
is as follows: normal order A{{a\ai}) and then replace the operators a\ and ctj 
by the fields one and a,i(t) respectively. This yields A(a(t)). We have adopted a 
path-integral formulation because of the large toolbox that goes with it. Besides, 
it formally provides us with a sort of dynamical partition function on which series 
expansion are straightforward. The efficiency of the mapping lies in the following 
observation: at e = 0, that is for independent particles, the Gaussian field theory 
exactly encodes Poissonian statistics. Hence the high-temperature (e — > 0) ex- 
pansion that will be performed next consists in expanding around the Poissonian 
distribution. For free particles, as expected, the stationary measure factorizes as a 
product of local Poissonian distributions, 

/-rii 

P(n)=U^fj (20) 

i 

Hence, for free particles, there is local equilibrium, and the density field = (rii) 
is identical to the fugacity field Q. This means that the properties of the system 
are the same as those that one would obtain in equilibrium if one imposed a space 
dependent fugacity Q. Here the explicit expression of the fugacity is given by, in 
the continuum limit, setting x — if L E [0,1], 

C(*) = Co + (Ci - Co)* + (7* - 0(1 - x)) + 0(L- 2 ) (21) 

where we have set Co = ^ and d — | as the fugacities of the reservoirs. The field 
theory described by So is free, which results in a(x, t) being a nonfluctuating field 
set to its average expression p(x) = ((x). This means in particular that, for e = 0, 

J VatVatAda^t)})^ 3 " = A({ Pl }) (22) 
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whatever the observable A. 



Another central quantity is the response function of the system to particle in- 
jection: the probability G(x, y;t — t') that there is a particle at x at time t given 
that there was one at y at time t' is given by 

G(x, y;r) = ^J2 *H^x) sin(n7n/)e^ 2n2r/L2 (23) 

The time integrated response function G(x, y) = J °° drG(x, y; r) will also be 
needed 

G(x, y) =L [x(l - y)Q(y - x) + y(l - x)G(x - y)} + 

l(l- x )(l- y )-l xy (24) 

- j^f (i* - Pi 1 - *)) - P0- - v)) + °( L ~ 2 ) 

In technical terms G is simply the free propagator of the theory. 



3.2 Effective free energy in the steady-state 

The probability P[{n;}] to observe a given occupation number configuration {n,} 
(or alternatively, a given proifle n(x)) is given by 

P[{n l }] = {\{5{n l -a\a i )) (25) 

i 

We further define the effective free energy J-\n\ of the profile n(x) by 

m = - lim (26) 

L^QO Li 

To access T we first pass to the generating function 

{n;} j 



then work directly on 



n[{*}] = - lim lnP p }] (28) 

L^OO Li 



which plays the role of an effective grand-potential. Our task is to compute 



P[{«}] = <e : 



Ei(*»-l)oi(*)\ 



(29) 



where the brackets denote the weighted path-integral defined by the action So- 
Using again a continuum limit notation, we find, as expected, that 



P[{z(x)}} 



£=0 



(30) 



from which we recover that the a steady-state probability distribution function is 
a product of local Poissonian distributions: 



rii 

i 

Going to a continuum notation we find 



71: 



(31) 



dx 



p(x) — n{x) + n(x) In 



n(x) 
p(x) 



(32) 



in agreement with the mathematically precise construction of the continuum limit 
0. At this stage, the present paragraph looks like a very technical rephrasing of 
simple properties. We are now ready, however, to attack the case of interacting 
systems. 



4 Zero Range Process 

4.1 Evolution operator and field-theory 

We will restrict our analysis to the case of pair repulsion, namely 

h(n) = en(n — 1) 
The master equation equation can be cast in the form 

<W)> 



dt 



(33) 



(34) 
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where for the rates specified in (1911 11121) . the evolution operator C £ takes the fol- 
lowing form 

-C £ =DJ2 E («! - «])e 2£a|ai a, 

i j nn of i 

+ a(l-aS) + <f(l-4) (35) 



7W 



l)e 2£a o ao a + /3(4 - l)e 2£a ^ aL a L 



One can of course recognize that at e = one recovers the free boson hamiltonian 
of the previous section. The corresponding action reads 



S £ [a, a] = — di(t) + J dr |^ ^dj(9 T aj + -D 



+ a(l - a ) + 7 («o - l)e^- 1)&oao ao 
+ 6(1 - a L ) + f3(a L - l)^ 6 "- 1 ^^ 



(36) 



Our calculations will be based on an expansion of S £ to first order in e; 



S £ = S + 2e / dt 



- a^a^a 2 + 7(a - l)a ao + P(a L - \)a L a\ 



+ 0(e 2 ) 
(37) 



4.2 Effective free energy, profile and correlations 

In order to evaluate 

P[z(x)] = (^/o^W*)- 1 )^.*)) (38) 

we rely on a cumulant expansion. This form is particularly well-suited for a cu- 
mulant expansion, which we write as 



P[^(rr)] = (exp 



1 r 1 n 

\{{^x 3 )~l)W {n \^-^n) + 

n>l n ' ^° j=l 

where is the n -point connected correlation function of the field a at equal 
times, in the steady state. In practice it is convenient to introduce the auxiliary 
fields (ft = a — 1 and </> = a — £. Note that 

p(x) = (n(x)} = (a(x)) = W {1 \x) (40) 



(39) 
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is so far undetermined. It will be useful to split into two pieces, denoted by 
Wfa! and Wtfz, corresponding to the local delta correlated piece in (namely 
the diagonal part) and the genuinely nonequilibrium contribution, respectively. In 
order to obtain one first directly computes the field connected correlation 
function to leading order in e: 

2e f°° f 1 / 

((f>{x 1 ,t 1 )(f)(x 2 ,t2)) c =— J dt J dy[C(y) 2 (dyG(x 1 ,y;t 1 -t)G(x 2 ,y;t 2 -t) 

+ G{x 1 ,y;t 1 -t)d 2 y G{x 2 ,y;t 2 -t))^ 

(41) 

where the usual contractions between a bared and a non-bared field were carried 
out. This yields, 

W®(x,y) = W®(x,y) = -2e( 2 (x)6{L{x - y)) (42) 

as if local equilibrium would hold (to this order in e at least). Note that the time 
depependent correlations are obtained by the same token. The profile for the ZRP 
reads 

p(x) = C(x) - 2e((x) 2 (43) 

Let us compare with what we would obtain in equilibrium for the distribution ©. 
The density-fugacity relationship and the local particle number fluctuations would 
read 

p = (n) = ( - 2e( 2 , (n 2 ) c = C-4< 2 (44) 

A quick glance at (1431421) shows that the zero range process, to first order in e and 
to leading order in L at least, appears to be consistent with local equilibrium. The 
existence of local equilibrium can actually be proved in general, independently of 
the explicit form of h(n) provided the transition rates are given by the expressions 
(1911 11121) . At the field theory level, it follows by inspection of the corresponding 
Feynman diagrams which indicates that local equilibrium holds to all orders in e 
due to the fact that the interaction arising from diffusion in the bulk is proportional 
to the Laplacian of the response field. Hence, by Wick's theorem, the interactions 
are proportional to the Laplacian of the propagator G which is a delta function in 
space. Thus spatial correlations cannot be built up and the steady state measure 
remains a product measure just as is the case for noninteracting particles d2T)l ). It 
should be noted, that the existence of local equilibrium can be more simply proved 
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JT51 by just substituting product measure into the steady state master equation and 
deriving an equation for the local fugacity (OTT) . 

For zero range processes the effective free energy can thus directly be ob- 
tained from the equilibrium distribution © in which one has substituted the local 
fugacity by its expression in terms of the local average density. 



5 Short range process 

5.1 Evolution operator and action 

As announced in 12.21 we now wish to study a different set of dynamical rules, 
those of the short-range process. For definiteness, we confine the present para- 
graph to the case of on-site pair repulsion, namely, h(n) — en(n — 1). One of the 
reasons for this choice is that, part from the 8 = limit which reduces to free par- 
ticles, the e — > oo limit exactly corresponds to the Symmetric Exclusion Process 
(SEP). Given that the scaling variable that appear in our expansions is ep, we may 
have the hope to connect our results with the low density behavior of the results 
obtained for the SEP. Again, for this dynamics, it is possible to write the master 
equation in the form of an imaginary time Schrodinger equation with an evolution 
operator 

i j nn of i 

+ a(l - al)e~ £ho + 5(1 - a[)e~ £hL (45) 
+ 7 (aS - l)e £ "°a + 0(4 - l)e £ ^a L 

In practice our analysis will be limited to the first nontrivial order in e, and the 
corresponding action reads 

S £ =So + e dt /](aj — a,i)(ai(ii — (Xjdj) 

- a(l - a )a a + 7(a - l)a a 2 Q (46) 

- 5(1 - a L )a L a L + (3(a L - l)a L a 2 L 

The expanded action (I46T) will be the starting point of our computations. 
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5.2 Effective free energy for onsite pair repulsion 

As explained earlier, we shall focus on the generating function of the probability 
distribution, 

P[z{x)\ = ( Q L foM^)-^,t)^ (47) 
Note that if one would start from a hamiltonian of the form 

h(n) = e g p n(n — l)...(n — p + 1) (48) 

p>2 

where the g p 's are order one constants, then the cumulant expansion would require 
to go as far as W^ p \ to leading order in e. Of course, at finite e, all cumulants are 
needed. In the present case of onsite pair repulsion, namely with h{n) = en{n—l) 
we are thus left with the computation of the first and second cumulant of a(x, t). 
We find that p(x) = (a(x, t)) has the following expression: 

p{x) =C(x) - 2e 7 (Co + C(0))C(0)G(x, 0) - 2e/?(Ci + C(l))C(l)G(x, 1) 

d y C{y)dyC{y)dyG{x,y) 

=(Co - 2eC 2 )(l -x) + (Ci - 2 £ Ci> + e(Ci - Co)Ml - x) 



+ e[Co(3 - 2x) + Ci(2 + 2x)]) 



(49) 



- P(l ~x)(l- e[Ci(l + 2x) + Co(4 - 2a;)]) j + 0{L- 2 ) 

It may be seen that 

p(x) = C(x) - 2e((x) 2 - e(Co - (i) 2 x(l - x) + 0(L- 1 ) (50) 

thus showing that there is no local equilibrium in the short range process. The 
density-density correlation function C^ 2 \x,y) = (n(x)n(y)) c takes the simple 
form 

CW(x,y) = (ax)-2eax) 2 )S(L(x-y)) + W^(x,y) (51) 

with 

WW(x,y)=-2e((x) 2 6(L(x-y)) 

- ^(Co - Ci) 2 [x(l - V)®(y -x)+y(l- x)Q(x - y)] (52) 
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(2) 

which now features a nonzero long range piece: 

w®(x, y) = - |(Co - CO 2 - y)®(y 1/(1 - x)e(x - v)\ (53) 

It is possible to return to the effective free energy 

T[n(x)} = l dx p(x) + ep(x) 2 — n(x) + n(x) In — — — - + en(x)(n(x) — 1) 

Jo L W + 2ep(a;) 2 



(54) 



The first brackets on the rhs of (154b corresponds to a system in local equilibrium 
with respect to an effective fugacity p+2ep 2 ; it already appeared for the zero range 
process. The integral in the second line illustrates the nonlocal long-range nature 
of the effective interactions in a nonequilibrium steady-state (W^ being negative, 
the corresponding contribution is positive, which expresses the repulsive nature of 
the effective interactions). The fluctuations of the total number of particles read 

AN 2 = AiY? 0C . eq . - L^( Po - Pl f (55) 

where AN 2 ^ eq denotes the fluctuations of a systems in local thermal equilib- 
rium with the same fugacity profile (such as the ZRP). The decrease of the global 
fluctuations is yet another consequence of the development of long-range anticor- 
relations. As a coincidence, note that setting in (1531) e = 1 yields exactly the 
symmetric exclusion process expression for this quantity Q. 



5.3 Effective free energy with onsite triplet repulsion 

Now we wish to explore the effect of varying the type of interaction by studying 
the case 

h(n) = en{n - l)(n - 2) (56) 

which assigns an energy cost to the site proportional to the number of triplets of 
particles are present at the site. We find that the profile is given by 

p(x) = ((x) - 3e((x) 3 -ex(l - a;) (Co - Ci) 2 (Co + Ci + Co(l - x) + &x) (57) 

V v ' 

local eq. 
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while the correlation function of the field reads 



W (2 \x, y) = - 6s((x) 3 5(L(x - y)) 



I2e 

- —(Co - Ci) 2 (Co + Ci) - y)e(y - x) + y{\ - x)Q(x - y)} 



(58) 



In d5Hl) the extra function w(x, y) has the explicit expression 



w(x,y) 
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7T 



4 ^— ' m? + n 2 

n,m>l 



(m + n) s 



n 



sin(n7rx) sin(m7n/) 



(59) 



where means a summation over ri and m of opposite parities. We have not 



n,m>l 



been able to come up with a closed expression for w (though one might exist). The 
function w vanishes at (1/2, 1/2) and does not contribute to the number of parti- 
cles fluctuations. Figure CO) shows a two-dimensional plot of w which reveals that 
w(x, y) conveys anticorrelations for x, y in the vicinity of the reservoir with the 
highest density, while positive correlations develop close to the reservoir impos- 
ing the lowest density. The function w stands for the deviation of the long-range 
component of the density correlation function with respect to the pair interaction 
case (whether at small e as in d52l . or at e — ► oo as computed by Spohn [ 16|). It 
is worth stressing that, to our knowledge, no microscopic model had, up to now, 
shown such strong deviations. 



The fluctuations of the total number of particles now read 



AN 2 = AN? oc . eq . 



2Le(Co-Ci) 2 (Co + Ci 



(60) 



where AA^ C eq denotes the fluctuations of a systems in local thermal equilibrium 
with the same fugacity profile (such as the ZRP). In the present case the third 
cumulant is of order e as well, and it has the formal expression and that 



W (3) {x 1 , x 2 , x 3 ) = -6eC(^i) 3 ^(^(^i - x 2 ))5(L(x 2 - x 3 )) 
- -J- / drd^(Co(l - z) + ( 1 z)G(xi,z;t)G(x 2 ,z;t)G(x 3 ,z;t) 



(61) 
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Figure 1: Plot of w(x, y) [see eq. (l59l) l as a function of x, y G [0, 1]. 



Denoting by the function appearing in the second line of the r.h.s. of Eq. doTl) 
allows to write the effective free energy in the form 



T[n(x)} = dx 



p + ep 3 — n + n In 



n 



p + 3ep 3 



dxdy 



nix) 



p{x) 



- 1 



wHH(x,y) 



n(y) 



+ en(n - l)(n - 2' 
n(y) _ 

p(y) 



p(y) 



- i 



n(z) 

W) 



1 I W^^j/^) 
(62) 



As one can see, the triplet repulsion generates not only two-body long ranged 
repulsive interactions, but also three body interactions, repulsive as well. This 
establishes a correspondence between the microscopies (the precise form of h{n)) 
and the effective interactions generated by the boundary drive. 



5.4 Competing interactions 

Suppose we now focus on a local interaction energy of the form 

h(n) = e [-n(n - 1) + Xn(n - l)(n - 2)] , A > 



(63) 
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where the first term on the r.h.s., which favors pair condensation, is competing 
with the second term on the r.h.s. which assigns a relative energy cost A to the 
piling up of three or more particles. To first order in e the density correlation 
function is merely a linear combination of that obtained for the pair and triplet 
interactions, namely 

C W( X <y ) = £(C - 0) 2 [(l - l2A(Co + Ci)Ml - V) + 12A(Co - (i)w(x, y) 

(64) 

We can clearly see that, at sufficiently low densities, positive correlations develop, 
thus providing a microscopic counterexample to the general belief that the driving 
of a current builds up anticorrelations (note, however, that Spohn [ 16 1 had left this 
possibility a priori open). 



6 Integrated current distribution 

Let Q(t) be the net number of particles which have jumped from the left reservoir 
into the system over the time interval [0, t] (counted positively for an actual jump 
from the reservoir into the system, and negatively for a jump out of the system 
to the reservoir). The physical motivations for studying the properties of Q are 
detailed by Lebowitz and Spohn [ 17 1 who showed that this quantity, defined for a 
Markov process (our boundary driven lattice gas with stochastic dynamics) plays 
a role analogous to the phase space contraction rate for dynamical systems (and 
for which Gallavotti and Cohen [ 18 1 proved their fluctuation theorem). Below, we 
shall concerned with the calculation of the distribution of Q(t). 



6.1 Free particles 

Following the procedure described in lElED, we construct a master equation for 
P{{ n i}i Qi t), the probability that the system is in state {rii] and with Q(t) = 
Q by explicitly separating those moves in phase space that increase Q by one, 
decrease it by one, or leave it unchanged. Introducing a state vector \ip(Q,t)} = 
^2r n .\ P({rii}, Q, t) | {rii}), we are left with an evolution equation of the form 

=~(Hi + H-i + H )MQ, t)) (65) 

dr 
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where the operators H±l increase/decrease Q by one, and where H leave Q 
unchanged. We are interested in the distribution function of Q(t) in the steady 
state, over large time intervals, 

p(Q,t) = {p\i/>(Q,t)) (66) 

or, more conveniently, its generating function, namely, 

+00 

p(z,t) = z Q P(Qit) (6V) 

Q=— 00 

It turns out that the generating function can conveniently be expressed in terms of 
a path-integral 

where, for noninteracting particles 

S 0>z [a, a) = - aj(t) + J dr [ ^ (^8^ + D ^(a^ ~ %) a ») 



+ a(l - za ) + 7(a - 2: 1 )a 
+ 5(1 - a L ) + (3(a L - l)a L 



(69) 



The action S 0j2 does not describe a stochastic process (unless z = 1 for which it 
reduces to So). However, it may be readily seen that 

p( Z ,t) = (e/o'dtW^-^-Ta-^^aoW)^ (70) 

where the brackets (..) z now denote an average with respect to the process gov- 
erned by So in which a is formally replaced with az. Using that ao(t), for free 
particles, is a nonfluctuating field with expression taken from (OTT) by changing a 
into az, 

a (t) = (oz - [1 - e(4( z + Ci)] (71) 

we find that, in the long time limit, 

\np(z,t) 1 z-1 
hm = fi(z) = — (C02; - Ci) (72) 

t^co t L z 
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It is a posteriori clear why subleading finite size corrections had to be kept all 
along. The function fj,{z) is the generating function for the cumulants of Q. It 
is worth commenting on two remarkable, yet expected, properties of /i(Co, Ci> z ). 
Namely, if we exchange the roles of the reservoirs, the distribution of Q is turned 
into that of — Q, hence 



MCo,Ci»«) = MCi,Co,2 l ) 



(73) 



But it also satisfies the Gallavotti-Cohen property (see ifTTl for a readable 
proof), namely 

/ i(Co,Ci^)=MCo,Ci,| L ) (74) 

which is best known when rephrased as follows. Let n(Q) be the large deviation 
function related to p(Q, t): 



lim 

t— »oo 



\np(Q,t) 



t 



(75) 



Besides, ir(Q) appears as the Legendre transform of fj,(z) with respect to In z: 

tt(Q) = max^{/i(Co, Ci, z) - In zQ} (76) 

hence 



li-niln^ 

t^OO t P[-Q,t) 



n(Q)-n(-Q)=\n^Q 
si 



(77) 



Of course this can a posteriori be verified on the explicit expression of ir(Q): 

CoCi 



<Q) =| + Vt + CoCi + 



-(Co + Ci)-Qln 



+ 



+ C0C1 



— + S0S1 + 2" 

Co 



(78) 



which is plotted on Fig. Q. The Gallavotti-Cohen theorem was shown lfT9lfT71l 
to hold under quite general conditions for nonequilibrium steady-states resulting 
from Markovian dynamics and we shall further comment upon it in subsection 
1631 
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Figure 2: Plot of ir(Q) in eq. (l78t for free particles with p = 1 ^> p\. Note the 
strong asymmetry of n. 



6.2 Short Range Process with pair repulsion 

We now repeat the strategy outlined in the previous paragraph, but for the interact- 
ing case, with the short range process dynamics. We split the evolution operator 
H e into three pieces, each describing moves that increase, decrease of leave Q 
unchanged. Then we pass to a path integral in terms of which we find 

«*.')= /aw* ™ 

Having found S £;Z we bring the calculation of p to that of a given observable 
with respect to the original stochastic process in which a is replaced with az: 

a(l — z)(l — ea a ) + 7(1 - z~ x )ao(l + ea Q a ) j) z 

(80) 

In order to evaluate the latter expectation value to leading order in e we rely on a 
cumulant expansion. Introducing the variable 

"(Co, Ci, z) = [(Co - eQ)z - (Ci - eg) - e(z - l)CoCi] (81) 

we find that 

L/jL — u- -uo 2 (82) 



p(z, t) = (exp 



dt 
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It is remarkable that the intermediate variable uo itself satisfies the two invariance 
properties that p is known to fulfill. Namely, the function uj(p , Pi, z) verifies the 
left-right symmetry 

w(Co,Ci > «)=w(Ci,Co,2" 1 ) (83) 
and the Gallavotti-Cohen property 

w(Co,Ci,*)="(Co,Ci,|^) C 84 ) 

That ui verifies (17317 41) implies directly that p possesses the same invariance prop- 
erties, as already identified in the SEP [6|. One invariance property that our result 
does not possess, in constrast to O, is the particle-hole symmetry. 

The effective small parameter of the expansion is ep, where p is the typical 
density, hence our small e expansion is seen to coincide with a low density expan- 
sion of the exact result obtained at e — > oo by [6|, 

Lp = lu - \uj 2 + ^-lu 3 + 0(u*) (85) 
6 45 

with 

Co Ci / ,\ Co c 

-z — (z — 1 



UJ 



i + Co~ i + Ci v ~ ^i + Coi + Cij (86) 

A natural question that arises next is how universal the result obtained in the exact 
calculation [6| is? How sensitive is it to varying the microscopic interactions. 
In order to sort out this issue we have performed a similar analysis for repulsive 
triplet interactions, with short range dynamics. 



6.3 Short Range Process with triplet repulsion 

In the same spirit as in the previous paragraph, we express the generating function 
p(z) as the expectation value of an exponential observable with respect to the 
lattice gas measure in which a is replaced with az. We rely again on a cumulant 
expansion, and after a rather tedious calculation, we arrive at the following results. 
We now define the auxiliary variable uj by 



w(Co,&>z) 



[(Co - e^)z - (Ci - eC!) - s(z - l)CoCi(Co + Ci 



x 



(<bs-Ci)(Co + Ci) 

2 z 



(87) 
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which also obeys (1731741 ) . We find that 



Lfj, — u — —lu 3 + 0{e 2 ) (88) 

This expression unambiguously points at a different distribution function for the 
integrated current. It further allows to connect the microscopies -the triplet repulsion- 
with the final form of fi. A p-body interaction would yield a first correction to \i 
of the form co p . Unfortunately we have not been able to come by a physical inter- 
pretation for the intermediate quantity uo. 

6.4 Zero Range Process with pair repulsion 

Finally we examine the case of the zero range process with pair repulsion, for 
which we know that the steady state distribution follows local equilibrium. It 
is then sufficient to start from the equilibrium expression for the free process in 
which one has subsituted the current with its appropriate expression. It is not hard 
to see, through a direct evaluation, that 

M = a _ 7 (n e 2£ ("°- 1 )) = (Co ~ 2<p 2 ) - (Ci ~ 2<i 2 ) „ Po - Pi (89) 
t L L 

which leads us to 

LMCo, Ci, z) = — [(Co - 2eCl)z - (d - 2eCD] (90) 

z 

This form is identical to that predicted by Bodineau and Derrida [12|) for zero 
range processes. The Gallavotti-Cohen relation, in our particular case, now takes 
the form 

MCo,Ci^) = MCo,Ci ^°~f^ -) (9D 

It is remarkable that at equilibrium, that is for p = p\, the current fluctuations 
(Q 2 ) is the same for the ZRP and the SRP, while higher order cumulants (that is 
the whole distribution) are different. 

6.5 Illustration of the additivity principle 

To conlude this section, we would like to illustrate how our explicit results fit 
into the general framework provided by Bodineau and Derrida |[T2l . In particular, 
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they have deduced, from a postulated additivity principle, a general expression 
for the integrated current distribution, from the sole knowledge of (Q) (p , pi) and 
(Q 2 ) c (po, Pi)- One of the consequences of their findings is a direct computation of 
the analog, in the Gallavotti-Cohen theorem, of the "entropy production rate" (as 
formally defined by Lebowitz and Spohn ifTTlO in these driven stochastic lattice 
gases. Defining, as in lfT2l. 



1 d 



t dp 

they have shown that 



(Q 2 )c(pip) 



D(P) = 7 — (Q)(Po,P) , *{P) = (92) 

P0=P I 



, : po D(p)i\ 

H(z)=n[2 &P-H--) (93) 



To leading order in e, it is straighforward to see from (1901) that, for the Zero Range 
Process, 

D(p) = 1, a(p) = 2p (94) 

For the short range process, 



for pair interaction: D(p) = 1 + 2ep, a(p) = 2p 
for triplet interaction: D(p) = 1 + Qep 2 , a(p) = 2p 



(95) 



thus leading, in both cases, to 



p(z) = Mtt) (96) 



Hence formula (193b is in perfect agreement with our explicit computations, not 
only for zero range processes, but also for the non trivial short range processes 
with pair or triplet repulsion. 



7 Final remarks 

We have shown on specific examples that driving a system out of equilibrium 
magnifies the differences in the underlying microscopic dynamics. Not only the 
density profiles are different, but some dynamical rules lead the system to a state 
of local equilibrium, while others let it develop long range effective interactions. 
In all cases, however, our explicit results for the integrated current distribution 
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provide support for the postulated additivity principle of Bodineau and Derrida 

G2 

We should emphasize that the above results have been worked out by setting 
up a formalism, based on path integrals, which provides an alternative both to 
exact solutions and to the fluctuating hydrodynamics approach. The path integral 
formalism fills a gap in the sense that it allows us to go directly from a microscopic 
formulation to macroscopic properties and, furthermore, it allows formulating ap- 
proximate approaches in nonequilibrium settings. Most of our results rely on a 
high temperature or virial like expansion, thus providing a intuitive parallel to the 
standard techniques of equilibrium statistical mechanics. 

We believe that several lines could now be explored. First, the present tool- 
box allows us to investigate the effect of additional space dimensions at little ex- 
tra formal cost (even though calculations will undoubtedly be more involved). 
This would help to clearly isolate those feature which are characteristic of one- 
dimensional systems from those that generalize to more realistic ones. An inter- 
esting problem arising in higher dimension is the interplay between longitudinal 
and transverse current fluctuations. Second, most of the quantities studied in the 
present work are time-independent, but the present toolbox makes possible the 
sudies of time-dependent quantities such as the time-dependent profile or effec- 
tive free energy. Third, a natural extension of our formalism would consist of 
establishing nonperturbative results in e. Much in the same way as in liquid state 
theory, infinite families of appropriately chosen Mayer diagrams can be summed 
up, and one could investigate which types of Feynman diagrams will contribute in 
building up the strongly nonlocal nature of the free energy functional found in JT). 
Finally, extending our results for more complex geometries with more than two 
particle reservoirs, as suggested in lfT2l or done in [20], should allow us to identify 
the universal emerging features and to bring the theory closer to experiments. 
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